Genome-Wide Characterization of Fennel (Anethum foeniculum) MiRNome and Identification of its Potential Targets in Homo sapiens and Arabidopsis thaliana: An Inter and Intra-species Computational Scrutiny.

MicroRNAs could be promising biomarkers for various diseases, and small RNA drugs have already been FDA approved for clinical use. This area of research is rapidly expanding and has significant potential for the future. Fennel (Anethum foeniculum) is a highly esteemed spice plant with economic and medicinal benefits, making it an invaluable asset in the pharmaceutical industry. To characterize the fennel miRNAs and their Arabidopsis thaliana and Homo sapience targets with functional enrichment analysis and human disease association. A homology-based computational approach characterized the MiRnome of the Anethum foeniculum genome and assessed its impact on Arabidopsis thaliana and Homo sapience transcriptomes. In addition, functional enrichment analysis was evaluated for both species' targets. Moreover, PPI network analysis, hub gene identification, and MD simulation analysis of the top hub node with fennel miRNA were incorporated. We have identified 100 miRNAs of fennel and their target genes, which include 2536 genes in Homo sapiens and 1314 genes in Arabidopsis thaliana. Functional enrichment analysis reveals 56 Arabidopsis thaliana targets of fennel miRNAs showed involvement in metabolic pathways. Highly enriched human KEGG pathways were associated with several diseases, especially cancer. The protein-protein interaction network of human targets determined the top ten nodes; from them, seven hub nodes, namely MAPK1, PIK3R1, STAT3, EGFR, KRAS, CDC42, and SMAD4, have shown their involvement in the pancreatic cancer pathway. Based on the Blast algorithm, 21 fennel miRNAs are homologs to 16 human miRNAs were predicted; from them, the CSPP1 target was a common target for afo-miR11117a-3p and has-miR-6880-5p homologs miRNAs. Our results are the first to report the 100 fennel miRNAs, and predictions for their endogenous and human target genes provide a basis for further understanding of Anethum foeniculum miRNAs and the biological processes and diseases with which they are associated.

Decoding the mechanism of andrographolide to combat hepatocellular carcinoma: a network pharmacology integrated molecular docking and dynamics approach.

HepatoCellular Carcinoma, being one of the most mortally convoluted malignancy with mounting number of occurrences across the world and being classified as the third most prevalent cause of cancer-associated mortalities and sixth most prevalent neoplasia. The active phytoconstituent andrographolide, derived from Andrographis paniculata is conveyed to reconcile a number of human ailments including various oncologies. However, the molecular mechanism underlying the anti-oncogenic effects of Andrographolide on HCC remains skeptical and unclear, emerging as a budding challenge for researchers and oncologists. The present study intends to analyze the underlying pharmacological mechanism of Andrographolide over HCC, established via assimilated approach of network pharmacology. Herein, the Network pharmacology stratagem was instigated to investigate potential HCC targets. The Andrographolide targets along with HCC targets were extracted from multiple databases. A total of 162 potential overlapping targets among HCC and Andrographolide were obtained and further subjected to gene ontology and Pathway enrichment analysis by employing OmicsBox and DAVID database, respectively. Subsequently, Protein-protein interaction network construction by Cytoscape software identified the top 10 hub nodes which were validated by survival and expression analysis. Further, the results derived from molecular docking and dynamic simulations by CB-Dock2 server and Desmond module (Schrodinger software) indicate ALB, CCND1, HIF1A, TNF, and VEGFA as potential Andrographolide related targets with high binding affinity and promising complex stability. Our findings not only reveal the antioncogenic role of andrographolide but also provide novel insights illuminating the identified targets as scientific foundation for anti-oncogenic clinical application of andrographolide in HCC therapeutics.Communicated by Ramaswamy H. Sarma.

MicroRNAs from Holarrhena pubescens stems: Identification by small RNA Sequencing and their Potential Contribution to Human Gene Targets.

Holarrhena pubescens is an effective medicinal plant from the Apocynaceae family, widely distributed over the Indian subcontinent and extensively used by Ayurveda and ethno-medicine systems without apparent side effects. We postulated that miRNAs, endogenous non-coding small RNAs that regulate gene expression at the post-transcriptional level, may, after ingestion into the human body, contribute to the medicinal properties of plants of this species by inducing regulated human gene expression to modulate. However, knowledge is scarce about miRNA in Holarrhena. In addition, to test the hypothesis on the potential pharmacological properties of miRNA, we performed a high-throughput sequencing analysis using the Next Generation Sequencing Illumina platform; 42,755,236 raw reads have been generated from H. pubescens stems from a library of small RNA isolated, identifying 687 known and 50 new miRNAs led. The novel H. pubescens miRNAs were predicted to regulate specific human genes, and subsequent annotations of gene functions suggested a possible role in various biological processes and signaling pathways, such as Wnt, MAPK, PI3K-Akt, and AMPK signaling pathways and endocytosis. The association of these putative targets with many diseases, including cancer, congenital malformations, nervous system disorders, and cystic fibrosis, has been demonstrated. The top hub proteins STAT3, MDM2, GSK3B, NANOG, IGF1, PRKCA, SNAP25, SRSF1, HTT, and SNCA show their interaction with human diseases, including cancer and cystic fibrosis. To our knowledge, this is the first report of uncovering H. pubescens miRNAs based on high-throughput sequencing and bioinformatics analysis. This study has provided new insight into a potential cross-species control of human gene expression. The potential for miRNA transfer should be evaluated as one possible mechanism of action to account for the beneficial properties of this valuable species.

Small RNA sequencing and identification of Andrographis paniculata miRNAs with potential cross­kingdom human gene targets.

Cross-species post-transcriptional regulatory potential of plant derived small non-coding microRNAs (miRNAs) has been well documented by plenteous studies. MicroRNAs are transferred to host cells via oral ingestion wherein they play a decisive role in regulation of host genes; thus, miRNAs have evolved as the nascent bioactive molecules imparting pharmacological values to traditionally used medicinal plants. The present study aims to investigate small RNA profiling in order to uncover the potential regulatory role of miRNAs derived from Andrographis paniculata, one of the most widely used herb by tribal communities for liver disorders and document the pharmacological properties of A. paniculata miRNAs. In this study, high-throughput sequencing method was used to generate raw data, ~ 60 million sequences were generated from A. paniculata leaves. Using computational tools and bioinformatics approach, analyses of 3,480,097 clean reads resulted in identification of 3440 known and 51 putative novel miRNAs regulating 1365 and 192 human genes respectively. Remarkably, the identified plausible novel miRNAs apa-miR-5, apa-miR-1, apa-miR-26, and apa-miR-30 are projected to target significant host genes including CDK6, IKBKB, TRAF3, CHD4, MECP2, and ADIPOQ. Subsequent annotations revealed probable involvement of the target genes in various pathways for instance p38-MAPK, AKT, AMPK, NF-Kß, ERK, WNT signalling, MYD88 dependant cascade, and pathways in cancer. Various diseases such as human papilloma virus infection, Alzheimer's, Non-alcoholic Fatty Liver, Alcoholic liver diseases, HepatoCellular Carcinoma (HCC), and numerous other cancers were predominantly found to be linked with target genes. Our findings postulate novel interpretations regarding modulation of human transcripts by A. paniculata miRNAs and exhibit the regulation of human diseases by plant-derived miRNAs. Though our study elucidates miRNAs as novel therapeutic agents, however, experimental validations for assessment of therapeutic potential of these miRNAs are still warranted.

MHC2AffyPred: A machine-learning approach to estimate affinity of MHC class II peptides based on structural interaction fingerprints.

Understanding how MHC class II (MHC-II) binding peptides with differing lengths exhibit specific interaction at the core and extended sites within the large MHC-II pocket is a very important aspect of immunological research for designing peptides. Certain efforts were made to generate peptide conformations amenable for MHC-II binding and calculate the binding energy of such complex formation but not directed toward developing a relationship between the peptide conformation in MHC-II structures and the binding affinity (BA) (IC50 ). We present here a machine-learning approach to calculate the BA of the peptides within the MHC-II pocket for HLA-DRA1, HLA-DRB1, HLA-DP, and HLA-DQ allotypes. Instead of generating ensembles of peptide conformations conventionally, the biased mode of conformations was created by considering the peptides in the crystal structures of pMHC-II complexes as the templates, followed by site-directed peptide docking. The structural interaction fingerprints generated from such docked pMHC-II structures along with the Moran autocorrelation descriptors were trained using a random forest regressor specific to each MHC-II peptide lengths (9-19). The entire workflow is automated using Linux shell and Perl scripts to promote the utilization of MHC2AffyPred program to any characterized MHC-II allotypes and is made for free access at https://github.com/SiddhiJani/MHC2AffyPred. The MHC2AffyPred attained better performance (correlation coefficient [CC] of .612-.898) than MHCII3D (.03-.594) and NetMHCIIpan-3.2 (.289-.692) programs in the HLA-DRA1, HLA-DRB1 types. Similarly, the MHC2AffyPred program achieved CC between .91 and .98 for HLA-DP and HLA-DQ peptides (13-mer to 17-mer). Further, a case study on MHC-II binding 15-mer peptides of severe acute respiratory syndrome coronavirus-2 showed very close competency in computing the IC50 values compared to the sequence-based NetMHCIIpan v3.2 and v4.0 programs with a correlation of .998 and .570, respectively.

A Comprehensive in vitro and in silico Assessment on Inhibition of CYP51B and Ergosterol Biosynthesis by Eugenol in Rhizopus oryzae.

Mucormycosis, also known as Zygomycosis, is a disease caused by invasive fungi, predominantly Rhizopus species belonging to the Order of Mucorales. Seeing from the chemistry perspective, heterocyclic compounds with an "azole" moiety are widely employed as antifungal agent for minimising the effect of mucormycosis as a prescribed treatment. These azoles serve as non-competitive inhibitors of fungal CYP51B by predominantly binding to its heme moiety, rendering its inhibition. However, long-term usage and abuse of azoles as antifungal medicines has resulted in drug resistance among certain fungal pathogens. Hence, there is an unmet need to find alternative therapeutic compounds. In present study, we used various in vitro tests to investigate the antifungal activity of eugenol against R. oryzae/R. arrhizus, including ergosterol quantification to test inhibition of ergosterol production mediated antifungal action. The minimum inhibitory concentration (MIC) value obtained for eugenol was 512 µg/ml with reduced ergosterol concentration of 77.11 ± 3.25% at MIC/2 concentration. Further, the molecular interactions of eugenol with fungal CYP51B were meticulously studied making use of proteomics in silico study including molecular docking and molecular dynamics simulations that showed eugenol to be strongly interacting with heme in an identical fashion to that shown by azole drugs (in this case, clotrimazole was evaluated). This is the first of a kind study showing the simulation study of eugenol with CYP51B of fungi. This inhibition results in ergosterol synthesis and is also studied and compared with keeping clotrimazole as a reference.

Unravelling the antifungal mode of action of curcumin by potential inhibition of CYP51B: A computational study validated in vitro on mucormycosis agent, Rhizopus oryzae.

Like human, fungi too are known to share lot of structural similarities amongst their CYPs (Cytochrome P450 super family of enzymes) which allows antifungal 'azole' compounds to interact with CYPs of human. Clotrimazole, an 'azole' antifungal drug, is a known inhibitor of fungal CYP named CYP51B. Curcumin, a phytochemical obtained from Curcuma longa has the ability to interact with several different human CYPs to induce inhibition. The sequence and the structural similarities amongst both human and fungal CYPs suggest a strong possibility for curcumin to interact with fungal CYP51B to behave like an antifungal agent. To test this hypothesis a study was designed involving mucormycosis agent, Rhizopus oryzae. The ability of curcumin to interact with fungal CYP51B was analysed computationally through molecular docking, MM-GBSA and Molecular Dynamics (MD) simulation assessment. Further, interaction profile for fungal CYP51B-curcumin was compared with human CYP3A4-curcumin, as there are published evidence describing curcumin as an inhibitor of human CYPs. Additionally, to validate in silico findings, an in vitro assay was performed to examine the antifungal potentials of curcumin on the R. oryzae. Conclusive results allow us to determine a plausible mode of action of curcumin to act as an antifungal against a mucormycosis agent.

Transcriptome-wide miRNA identification of Bacopa monnieri: a cross-kingdom approach.

Bacopa monnieri known as 'Brahmi' is a well-known medicinal plant belonging to Scrophulariaceae family for its nootropic properties. To the best of our knowledge, no characterization data is available on the potential role of micro RNAs (miRNAs) from this plant till date. We present here the first report of computational characterizations of miRNAs from B. monnieri. Owing to the high conservation of miRNAs in nature, new and potential miRNAs can be identified in plants using in silico techniques. Using the plant miRNA sequences present in the miRBase repository, a total of 12 miRNAs were identified from B. monnieri which pertained to 11 miRNA families from the shoot and root transcriptome data. Furthermore, gene ontology analysis of the identified 68 human target genes exhibited significance in various biological processes. These human target genes were associated with signaling pathways like NF-kB and MAPK with TRAF2, CBX1, IL1B, ITGA4 and ITGB1BP1 as the top five hub nodes. This cross-kingdom study provides initial insights about the potential of miRNA-mediated cross-kingdom regulation and unravels the essential target genes of human with implications in numerous human diseases including cancer.

Molecular interaction of selected phytochemicals under the charged environment of Plasmodium falciparum chloroquine resistance transporter (PfCRT) model.

Phytochemicals of Catharanthus roseus Linn. and Tylophora indica have been known for their inhibition of malarial parasite, Plasmodium falciparum in cell culture. Resistance to chloroquine (CQ), a widely used antimalarial drug, is due to the CQ resistance transporter (CRT) system. The present study deals with computational modeling of Plasmodium falciparum chloroquine resistance transporter (PfCRT) protein and development of charged environment to mimic a condition of resistance. The model of PfCRT was developed using Protein homology/analogy engine (PHYRE ver 0.2) and was validated based on the results obtained using PSI-PRED. Subsequently, molecular interactions of selected phytochemicals extracted from C. roseus Linn. and T. indica were studied using multiple-iterated genetic algorithm-based docking protocol in order to investigate the translocation of these legends across the PfCRT protein. Further, molecular dynamics studies exhibiting interaction energy estimates of these compounds within the active site of the protein showed that compounds are more selective toward PfCRT. Clusters of conformations with the free energy of binding were estimated which clearly demonstrated the potential channel and by this means the translocation across the PfCRT is anticipated.

Structural insights into the theoretical model of Plasmodium falciparum NADH dehydrogenase and its interaction with artemisinin and derivatives: towards global health therapeutics.

It is a continuing quest to uncover the principal molecular targets of malarial parasites to understand the antimalarial activity and mechanism of action of artemisinin, a potent antimalarial. A series of parasite proteins are experimentally validated as potential targets, such as translationally controlled tumor protein (TCTP) and sarco/endoplasmic reticulum membrane calcium ATP-ase (SERCA). The present study addressed the development of a theoretical model of Plasmodium falciparum NADH dehydrogenase with inference from artemisinin in vivo inhibitory activity. We report here the predicted binding modes of artemisinin and its derivatives. The modeled protein resembled the structural architecture of flavoproteins and oxidoreductases, consisting of two Rossmann folds and dedicated binding sites for its cofactors. Docked poses of the ligand dataset revealed its interactions at or near the si face, indicating being activated. This may aid in generation of reactive oxygen species, thereby disrupting the membrane potential of parasite mitochondria and leading to the clearance from the blood. These observations open up new strategies for development of novel therapeutics, or improvement of existing pharmacotherapies against malaria, a major burden for global health.

A Computational Approach towards the Understanding of Plasmodium falciparum Multidrug Resistance Protein 1.

The emergence of drug resistance in Plasmodium falciparum tremendously affected the chemotherapy worldwide while the intense distribution of chloroquine-resistant strains in most of the endemic areas added more complications in the treatment of malaria. The situation has even worsened by the lack of molecular mechanism to understand the resistance conferred by Plasmodia species. Recent studies have suggested the association of antimalarial resistance with P. falciparum multidrug resistance protein 1 (PfMDR1), an ATP-binding cassette (ABC) transporter and a homologue of human P-glycoprotein 1 (P-gp1). The present study deals about the development of PfMDR1 computational model and the model of substrate transport across PfMDR1 with insights derived from conformations relative to inward- and outward-facing topologies that switch on/off the transportation system. Comparison of ATP docked positions and its structural motif binding properties were found to be similar among other ATPases, and thereby contributes to NBD domains dimerization, a unique structural agreement noticed in Mus musculus Pgp and Escherichia coli MDR transporter homolog (MsbA). The interaction of leading antimalarials and phytochemicals within the active pocket of both wild-type and mutant-type PfMDR1 demonstrated the mode of binding and provided insights of less binding affinity thereby contributing to parasite's resistance mechanism.

Evolutionary and molecular aspects of Indian tomato leaf curl virus coat protein.

Tomato leaf curl disease (ToLCD) is manifested by yellowing of leaf lamina with upward leaf curl, leaf distortion, shrinking of the leaf surface, and stunted plant growth caused by tomato leaf curl virus (ToLCV). In the present study, using computational methods we explored the evolutionary and molecular prospects of viral coat protein derived from an isolate of Vadodara district, Gujarat (ToLCGV-[Vad]), India. We found that the amino acids in coat protein required for systemic infection, viral particle formation, and insect transmission to host cells were conserved amongst Indian strains. Phylogenetic studies on Indian ToLCV coat proteins showed evolutionary compatibility with other viral taxa. Modeling of coat protein revealed a topology similar to characteristic Geminate viral particle consisting of antiparallel ß-barrel motif with N-terminus α-helix. The molecular interaction of coat protein with the viral DNA required for encapsidation and nuclear shuttling was investigated through sequence- and structure-based approaches. We further emphasized the role of loops in coat protein structure as molecular recognition interface.

Epitope-based immunoinformatics and molecular docking studies of nucleocapsid protein and ovarian tumor domain of crimean-congo hemorrhagic Fever virus.

Crimean-Congo hemorrhagic fever virus (CCHFV), the fatal human pathogen is transmitted to humans by tick bite, or exposure to infected blood or tissues of infected livestock. The CCHFV genome consists of three RNA segments namely, S, M, and L. The unusual large viral L protein has an ovarian tumor (OTU) protease domain located in the N terminus. It is likely that the protein may be autoproteolytically cleaved to generate the active virus L polymerase with additional functions. Identification of the epitope regions of the virus is important for the diagnosis, phylogeny studies, and drug discovery. Early diagnosis and treatment of CCHF infection is critical to the survival of patients and the control of the disease. In this study, we undertook different in silico approaches using molecular docking and immunoinformatics tools to predict epitopes which can be helpful for vaccine designing. Small molecule ligands against OTU domain and protein-protein interaction between a viral and a host protein have been studied using docking tools.